SYSTEM AND METHOD FOR GEOREFERENCING DIGITAL 
RASTER MAPS 

CROSS REFERENCE TO RELATED APPLICATIONS: 

This application is a continuation of and claims priority from pending 

5 U.S. Patent Application "System and Method for Georeferencing Digital Raster 

Maps" (09/537,849), filed March 29, 2000. Furthermore, this apphcation is 

related to and claims priority from pending applications: "System and Method 

for Performing Flood Zone Certifications" (09/537,161), filed March 29, 2000 

and "System and Method for Synchronizing Raster And Vector Map Images" 

10 (09/537,162), filed March 29, 2000 which are hereby incorporated by reference. 

BACKGROUND OF THE INVENTION 
Technical Field: 

The present invention generally relates to geographic map processing and 
in particular to a system and method for associating digital raster maps with 

1 5 known geographic features. Still more particularly, the present invention relates 
to a system and method for associating specific points on digital raster maps with 
a geographic coordinate system. 
Description of the Related Art: 

A digital raster map is a computerized map image that resembles a 

20 conventional paper map in that it presents an image of the mapped area, but has 

no additional underlying data associated with the features of the map. A raster 
map is typically created by scanning a conventional paper map, and is a grid- 
based map composed of pixels (or dots) of color or black & white. Each pixel in 
the grid has can be referenced by the pixel coordinates, and has only one value, 

25 which indicates the color of that pixel. Raster images are commonly referred to 

as "bit mapped." 

A vector map uses lines and polygons, rather than pixels, to describe an 
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image. Unlike a raster map, in which the map image is simply stored as a grid of 
pixels, when a vector map is displayed, it is drawn from a set of underlying data. 

The vector map is created through the placements of nodes on a plane and 
connecting those nodes with lines. Vector lines can be attributed with tables of 
5 data such as elevations, values, names or other information relative to the line. 

Vector data can be displayed in three dimensions if the lines are attributed with 
z values, modified or changed relative to user need, or layered to allow for turning 
off and on the viewing of different information. 

Because of their feature attribution properties, they are particularly useful 
10 for displaying geographic data. Vector maps are used to display boundaries or 

lines that denote the position and extent of features, such as county boundaries or 
lines denoting stream and river systems. It is also very easy to view or manipulate 
the data underlying a vector map, for example to view or change the elevation of 
a feature. 

15 Also because vector maps are commonly drawn from a geographic data 

set describing the area shown, they are very easily, and generally inherently, 
georeferenced. Georeferencing is the process of relating source coordinates to 
referenced geographic coordinates., which are typically in standard 
latitude/longitude. An image or a vector file is georeferenced to be used within 

20 a mapping/geographic environment. In a vector map, the data from which the 

map is drawn will typically already include a geographic coordinate set. 

Modem GIS systems normally make use of digital vector based map 
information. However, a vast legacy of paper-based map information exists. It 
is very expensive and time consuming to convert all of the information on these 
25 paper maps over to a digital vector format. In many cases the scope and expense 

of such conversions render them completely impractical. However, even when 
a complete conversion to digital vector-based format is not possible, it is still 
possible to obtain some of the benefits of computerized map systems, first by 
converting the paper maps to digital raster maps (by scanning them), and then by 
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georeferencing the raster image. After georeferencing, there should be a clear 
relationship between the pixel coordinates in the raster map, and the geographic 
coordinates of the feature represented by that pixel. It would be desirable to 
provide a system and method for georeferencing a raster map by associating 
5 points on that map with corresponding points on a previously-georeferenced 

vector map, or with reference points designated in the raster map which have 
known latitude and longitude. 
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SUMMARY OF THE INVENTION 

It is therefore one object of the present invention to provide improved 
geographic map processing. 

It is another object of the present invention to provide an improved 
5 system and method for associating digital raster maps with known geographic 

features. 

It is yet another object of the present invention to provide an improved 
system and method for associating specific points on digital raster maps with a 
geographic coordinate system. 

10 The foregoing objects are achieved as is now described. The preferred 

embodiment provides a system and method for georeferencing digital raster maps 
by associating specific points on the raster map with corresponding points on a 
georeferenced vector map or another georeferenced raster map. According to the 
preferred embodiment, a raster map and a corresponding vector map are 

1 5 simultaneously displayed to a user. The user then locates a common geographic 

point or feature on each map, and marks each of them as a unique point-pair. 
When the user has marked at least two point-pairs, the system then computes a 
georeferencing function, based on the pixel-coordinates of the points marked on 
the raster map and the corresponding geographic coordinates of the points on the 

20 vector map. Thereafter the geographic coordinates of any point on the map may 

be easily computed. The preferred embodiment provides that as more point-pairs 
are identified, the georeferencing function is modified for increased accuracy. 

The above as well as additional objectives, features, and advantages of the 
present invention will become apparent in the following detailed written 
25 description. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

The novel features believed characteristic of the invention are set forth in 
the appended claims. The invention itself however, as well as a preferred mode 
of use, further objects and advantages thereof, will best be understood by 
reference to the following detailed description of an illustrative embodiment 
when read in conjunction with the accompanying drawings, wherein: 

Figure 1 depicts a data processing system in accordance with a preferred 
embodiment of the present invention; 

Figure 2 is an image of an exemplary raster map, in accordance with the 
preferred embodiment; 

Figure 3 is an image of an exemplary vector map, corresponding to the 
raster map of Figure 2, in accordance with a preferred embodiment of the present 
invention; and 

Figure 4 is a flowchart of a process in accordance with a preferred 
embodiment of the present invention. 



108344.00011: 2604976 



DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 



With reference now to the figures, and in particular with reference to 
Figure 1, a block diagram of a data processing system in which a preferred 
embodiment of the present invention may be implemented is depicted. Data 
5 processing system 100 includes processor 102 and associated L2 Cache 104, 

which in the exemplary embodiment is connected in tvtm to a system bus 106. 
System memory 108 is connected to system bus 106, and maybe read firom and 
written to by processor 102. 

Also connected to system bus 106 is I/O bus bridge 110. hi the exemplary 

1 0 embodiment, data processing system 100 includes graphics adapter 118 connected 

to bus 106, receiving user interface information for display 120. Peripheral 
devices such as nonvolatile storage 114, which may be a hard disk drive, and 
keyboard/pointing device 116, which may include a conventional mouse, a 
trackball, or the like, are connected to I/O bus 112. 

1 5 The exemplary embodiment shown in Figure 1 is provided solely for the 

purposes of explaining the invention and those skilled in the art will recognize 
that numerous variations are possible, both in form and function. For instance, 
data processing system 100 might also include a compact disk read only memory 
(CD-ROM) or digital video disk (DVD) drive, a sound card and audio speakers, 

20 and numerous other optional components. All such variations are believed to be 

within the spirit and scope of the present invention. Data processing system 100 
is provided solely as an example for the purposes of explanation and is not 
intended to imply architectural limitations. 

The preferred embodiment provides a system and method for 

25 georeferencing digital raster maps by associating specific points on the raster map 

with corresponding points on a georeferenced vector map or another 
georeferenced raster map. According to the preferred embodiment, a raster map 
and a corresponding vector map are simultaneously displayed to a user. The user 
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then locates a common geographic point or feature on each map, and marks each 
of them as a xonique point pair. When the user has marked at least two point-pairs, 
the system then computes a georeferencing function, based on the pixel- 
coordinates of the points marked on the raster map and the corresponding 
5 geographic coordinates of tiie points on the vector map. Thereafter the geographic 

coordinates of any point on the map may be easily computed. The preferred 
embodiment provides that as more point-pairs are identified, the georeferencing 
function is modified for increased accuracy. 

Figure 2 is an exemplary raster map, in accordance with the preferred 

10 embodiment. This exemplary map shows a scanned image firom a Federal 

Emergency Management Agency (FEMA) paper map. This raster image shows 
a land area with flood zone indications, but would, in a computer system, contain 
no underlying data regarding the area shown. 

Figure 3 is an exemplary vector map, corresponding to the raster map of 

1 5 Figure 2, in accordance with a preferred embodiment of the present invention. 

This map shows the same area as the map in Figure 2, but is created by a 
computer system fi-om a database describing the locations of features such as the 
streets shown. Typically, each feature shown on a vector map such as this will 
already be georeferenced, in that the geographic coordinates of each feature will 

20 also be recorded in the underlying data 

A digital map image is considered georeferenced if a pair of mathematical 
functions, f, and g, have been defined that can be used to convert back and forth 
between the coordinates of the map image (as defined by the pixels of the image) 
and the corresponding longitude and latitude of the location of that point. That is, 

25 f and g do the following: 

1 . If (x, y) represents a location on the digital map image, then f(x, 
y) HJjon, Lat) represents the longitude and latitude of the corresponding physical 
location. 
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2. If (Lon, Lat) represents a physical location that lies within the 
region covered by the map, then g (Lon, Lat) = (x, y) represents the point on the 
digital map image that corresponds to that longitude and latitude. 

Here, x and y represent the natural internal coordinate system of the map 
5 image. Typically, as described above, a digital raster map image uses the pixels 

of its image as a natural coordinate matrix. However, in most cases, a vector- 
based map image uses longitude and latitude as its intemal coordinate system; if 
so, it can be considered to be trivially georeferenced already. Therefore, the 
functions f() and g ( ) above are non-trivial georeferencing functions 

1 0 required to convert back and forth between coordinate systems. 

Following is a description of a georeferencing process in accordance with 
the preferred embodiment: 

First, the raster map to be georeferenced (referred to as Mapl) is 
displayed on a computer monitor in conjunction with another previously 

1 5 georeferenced map (referred to as Map2) (step 400). Typically, Map2 will be a 

vector map, but could also be another raster map, for which georeferencing 
functions are already known. Mapl and Map2 can be independently manipulated 
for rotation, position and scale. The user will manipulate the maps until the region 
shown in Mapl is contained within the region shown by Map2. 

20 Next, pairs of points representing identical geographic locations on Mapl 

and Map2 are identified, and the user uses a keyboard, mouse, or other pointing 
device to mark these point-pairs (step 405). For each point-pair, one point is 
marked on Mapl and the other point is marked at the corresponding location on 
Map2. Examples of common georeferencing point-pairs mclude street 

25 intersections, places where rivers cross roads, mountain peaks, building locations, 

and other readily identifiable map features. If any points on Mapl have points of 
known longitude and latitude already marked on the map, then a georeferencing 
point-pair can be obtained by marking that point on Mapl and moving to the point 
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with identical latitude and longitude on Map2 and marking it as the other half of 
the point-pair. Alternatively, once the known points are marked with the known 
latitude and longitude, they can be used to determine the georeferencing functions 
without requiring a corresponding mark on Map2. 

5 After a point-pair is marked, the x-y point on Mapl is assigned the 

corresponding latitude and longitude of its matching point on Map2 (step 410). 

When two georeferencing point-pairs are marked, the system will 
compute georeference functions for Mapl based on a linear transformation that 
allows an arbitrary rotation and predefined scaling, between the x-y domain and 

1 0 the latitude- longitude domain (step 415). 

To improve accuracy and allow increased generaUty, more georeferencing 
point-pairs will frequently be desired. To facilitate the process of gathering these 
additional georeferencing point-pairs, Mapl and Map2 are now selectively 
synchronized. This means that as either map is zoomed, panned, scrolled, or 

1 5 otherwise caused to display a different region, then the other map automatically 

does the same. This greatly increases the ease with which georeferencing point- 
pairs can be identified and marked on the maps, since when the user of the system 
locates an identifiable feature on one map, the other map will automatically be 
showing very nearly the same location on the other map. 

20 Thereafter, when the user marks a georeferencing point on one of the 

maps, then the system automatically marks the corresponding point on the other 
map, based on its current georeferencing function (step 420). The user can then 
make corrections to the location of a point marked by the system in this way (step 
425). The user may enter as many point-pairs as are desired; each additional 

25 point-pair will increase the accuracy of the georeferencing process. 

When three or more georeferencing point-pairs are determined, the system 
computes georeference fimctions for Mapl based on a completely general linear 
transformation between the x-y domain and the latitude-longitude domain. 
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When four or more georeferencing point-pairs are determined, the general 
linear georeferencing functions are over-determined. This means that more than 
the required amount of information to compute the general linear georeferencing 
functions is available, but that it is not, in general, completely consistent. The 
5 system use the extra information contained in the additional georeferencing points 

to provide validation checks to protect against the possibihty that some of the data 
points may be inaccurate (step 430). Points that deviate excessively with respect 
to a calculated standard error are presumed to be inaccurate and are omitted from 
the calculation of the georeferencing functions. Note that as new point are added, 
1 0 the system also rechecks points previously marked as inconsistent, to determine 

if those points should now be considered when recomputing the georeferencing 
functions. 

The additional point-pairs allow the system to compute the general linear 
georeferencing functions which best fit the combination of all the available 
1 5 georeferencmg point-pairs, so that the active georeferencing functions are revised 

with each new point-pair (step 435). According to the preferred embodiment, this 
is done by using a "least square" parameter fitting operation. 

The user may then proceed o enter the next point-pair (step 440) . When 
the user is finished, the system stores the active georeferencing functions with the 
20 raster-map (step is 445). At this time, the raster map is considered fully 

georeferenced. When accessed at any future time, the system may simply retrieve 
the georeferencing functions, and apply them to find the latitude and longitude of 
any point on the raster map. 

The process of determining he georeferencing function set fi:om a set of 
25 point-pairs is beheved to be within the ability of one of ordinary skill in the art. 

The specific ^proach used by the system an method of the preferred embodunent 
is discussed below. 

At any given point in the georeferencing process a set of points P , has 
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been specified. Let / e P denote one of the points of P . Associated value x, , y,, 
Lat, and Lon, refer to thejc and y coordinates of the point on the bitmap, and 
the latitude and longitude of the point on the earth's surface. We also define 
A c P to be the subset of "active" points, i.e. points which are currently being 
5 used in the calculation of the georeferencing fiuictions. Assume that we have 

selected a parametric family of fiinctions F. Functions in this family map fi-om 
{x, y) to {ion, Lat) . Any fimction / e calF has components /j {x, y) = Lon, and 
fi y) = Lat . From within this family we seek the specific fimction / , which 
(in a sense to be defined below) comes closest to making the following system of 
10 equations true: 

/(x, ,y,) = {Lon, , Lat, ) for i e A 

Once determined, / will b the georeferencing fimction which is used to 
compute corresponding latitude and longitude values, {Lon, Lat) for any point, 
{x,y) on the bitmap. There are any number of possible ways to define the 
1 5 fimction that "comes closest to making (1) true." We shall follow a "least squares" 

approach also known in mathematics as an approach. This approach seeks to 
find the fimction, / , which minimizes the sum of the squared differences 
between the actual and the predicted values of latitude and longitude, hi other 
words, fi-om among all the fimctions / e F , / is the one which minimizes: 

20 SSE = Y!if,{x„yyLon,Y +{f^{x„y,)- Lat,Y\ (2) 

Among various alternative methods for choosing the fimction /are 
choosing it so t ht it minimizes the simi of absolute errors (rather than squared 
errors, or so that it minimizes the largest error. Other criteria are also possible. 
In what follows we shall describe methods for computing georeferencing 
25 fimctions based on both a "rotational linear" parametric family of fimctions and 
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a "general linear" family of functions. A knowledgeable practitioner will readily 
perceive that this approach may be easily generalized to other parametric famihes 
of functions. Obvious candidates include nonlinear parametric families of 
functions arising from the projection of the earth's surface onto a flat map. 

5 The General Linear Case: Li the general linear case, we let F be the set 

of all possible hnear transformations which map from (x, y) to {Lon, Lat) . Thus 

for some choice of parameters a^^ , a^2> ^21' ^22* ^i' ^2- If the 
region covered by the map to be georeferenced is not too large, then this family 

10 of functions will contain a suitable function , / whose total error is quite small. 

In the case where the map to be georeferenced covers a larger area than this, then 
the curvature of the earth must be taken into account and F is not a suitable 
family of functions, hi such a case, nonlinear functions must be used as 
mentioned above. We shall not pursue that case further, since it is a 

1 5 straightforward extension of the procedures used in the linear case. 

To find / we seek the parameters which minimize 
SSE = 2] {a^^Xi + «i2>', +b,- Lon, f + (gji^^, + ^22)^, + 63 " f The 

parameter values which minimize this expression are foimd by solving the 
following two independent systems of linear equations: 

20 

n 
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^JC; Lon, 
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Y^yiLon, 



(5a) 
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S^/ Z^' E^' 



i€A 



(5b) 



These systems can be easily solve by well known methods, such as 
Gaussian Elimination, or LU factorization. The solutions yield the desired values 
of a„ , a,2 , ' ^22' 4 ' 3nd ftj • It should be noted that equations (5a) and (5b) 
do not have a unique solution unless three or more non-colinear points are 
contained in A. Generally speaking, then, it requires 3 points to choose a 
georeferencing function from the family of general linear transformations. When 
there are four points or more, it is possible to compute a standard deviation of 
errors using the formula: 



^ (qijJC, + 5i2 J, + 6i - Lon, y + {a2iX, +8^- Lat, J 



(6) 



is an estimator for the amount of error to be expected between actual 
1 5 and predicted latitude and longitude values. 

Note that the inverse georeferencing function, that maps from (Lon, 
Lat) to (x,y) is readily obtained now, by inverting the function / . Having done 
this, it is possible to compute a similar standard deviation of error for which 
is an estimate of error as measured in the bitmap coordinates. 
20 The Rotational Linear Case: In the, so called, rotational linear case, we 

let Fbe the set of all possible linear transformations which map from (x,>')to 
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{Lon,Lat) and which also allow for a known scale change caused by the 
difference in distance per degree of latitude and per degree of longitude at any 
point not on the equator. Thus 



-A' 



(7) 



for some 

choice of parameters , A» A» where / =cos(la/) is 
computed firom the estimated latitude of the map to be georeferenced. 
Clearly, this is a special case of the more general linear transform 
described above, where we take 

In theory, if the region covered by the map to be georeferenced is not too 
large, then this family of functions will contain a suitable function, / 
whose total error is quite small. In the case where the map to be 
georeferenced covers a larger area than this, then the curvature of the 
earth must be taken into account and F is not a suitable family of 
functions. In such a case, nonlinear functions must be used as noted 
above. In practice, the general linear transform will yield results which are 
somewhat superior to the rotational linear transform due to the former's 
innate ability to use arbitrary scaling factors. These arbitrary scaling 
factors might compensate, for example, for scaling flaws in the digital 
scanner used to scan in the maps. The primary advantage of the rotational 
linear transform is that it can be computed with only two reference points 
(unlike the general linear transform which requires three reference points). 
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To find / we seek, as before, the parameters which minimize SSE as 
defined in equation (4). 

The parameter values which minimize this expression are found by 
solving the following system of linear equations: 
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These systems can easily solved by well known methods, such as 
Gaussian Elimination, or LU factorization. The solutions yield the desire values 
10 of , , >§i , and y^j which in turn yield the desired values for , > ^21 ' 

It should be noted that equation (8) does not have a unique solution unless 
two or m re points are contained in A. Generally speaking, then it requires two 
1 5 points to determine a georeferencing function fi-om the family of rotational linear 

transformations. When there are three points or more, it is possible to compute 
a standard deviation of error, s usmg the formula: 
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s is an estimator for the amount of error to be expected between actual 
and predicted latitude and longitude values. 

Note that the inverse georeferencing function /"^ that maps from (Lon, 
Lat) to (jc,>^) is readily obtained now, by inverting the function / . Having done 
this, it is possible to compute a similar standard deviation of error for which 
is an estimate of error as measured in the bitmap coordinates. 

Automatic Error Detection and Handling When individual points are 
being assigned x, y, Lon, and Lat values, there is always a potential for error. To 
reduce the risk of incorrect georeferencing resulting from such errors, certain 
error handling procedures are built into the georeferencing process. The 
fundamental concept is that of detecting a "bad" point and then removing it from 
the set of active points, A Note that removing point from A will not delete the 
information associate, with that point, but it will cause the georeferencing 
parameters to be completely uninfluenced by that point. We o not wish to remove 
the point entirely, since it may be determined at a later stage of the 
georeferencing, that the point was not really bad at all, and should be-used in the 
georeferencing calculation. This will be clarified shortly. 



20 Detecting Bad Points The following steps outline the bad point detection 

process using the general linear transform approach to georeferencing. 

1 . Begin by placing all existing points into the active set, A 

2. If there are fewer than five active points then you are done. 
Otherwise, for each of the currently active points in turn, move it (call 

25 it point k for the sake of convenience) temporarily out of the active set, and 
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then calculate the resulting inverse georeferencing function (call it g^*^) and 
its corresponding SSE^ . Also, calculate the difference betweai the predicted 
value and the actual value 5/, = ^g^^'^Lon/^ , Lat^. ) - (x^ , )| . Make a note of the 

values, Sj^and ■ Return point k to the active set (and move on to the next 

value of k . 

3. From among the results found in step 2 above, find the point 
k, with the largest value of ^y^^^ )Cyand5^ which also satisfies 

^y^l^^ }c^andSi^)c 2, where Cy, andcj are some constants which are set 

according to the general level of accuracy to be expected on the particular type of 
map which is being georeferenced, the current number of active points, and the 
dots per inch of the scanned image. If there is such a point then mark it as bad (by 
removing it firom the active set) and return to step 2 above. Otherwise you are 
done. 

There are several things to note about this procedure. One is that it 
allowing the value of and to change with the number of active points, makes 
it possible for the georeferencing system and method to utilize points which it 
might originally determine bad or inconsistent after a large enough sample of 
points has been gathered to make it clear that a lesser level of accuracy is all that 
can be achieved on this map. Another observation is that using this procedure it 
is impossible to reduce the number of active points down to less than four (unless 
you started with less than 4 in which case this procedure does not apply at all). 
This scheme means that as each new point is added, all points determined so far 
are considered, even those which had previously been marked bad. Thus early 
"misjudgments" on the part of the system can be corrected later, in Ught of new 
point information. 

The same bad point detection process, can also be implemented using the 

108344.00011: 2604976 J 7 



rotational linear transform approach. In this case the method is capable of 
reducing the number of active points down as low as three (rather than four for 
the general linear transform approach outlined above). This can be useful when 
dealing with small sets of active points. 

When the system is georeferencing raster map images that cover a large 
enough area so that the (nonlinear) curvature of the earth is a source of significant 
error, the system can calculate nonlinear georeferencing functions whose form 
corresponds to the map projection that was used to create Mapl. 

A specific example of the operation and application of the preferred 
georeferencing method may be shown with reference to the "Flood Zone 
Determination" business. The Federal Emergency Management Agency (FEMA) 
publishes a library of tens of thousands of paper maps showing various types of 
flood zones and their locations in the United States. A flood zone determination 
on a property is fi-equently done in the following way: 

1 . The address of the property is examined, and the 
location of the property is determined (perhaps through the 
use of a geocoding system, orb examining an available 
street map). 

2. A map analyst attempts to determine which of the 
many thousands of FEMA flood maps will contain this 
property. 

3. The map analyst goes to a map storage area and 
retrieves the desired map, often examining several maps 
before making a final selection. 

4. Having retrieved the paper map, the map analyst 
next determines where, precisely, the property is located on 
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the map. 

5. Finally, the map analyst examines flood zone 
notations on the map at the property's location in order to 
determine its flood- zone status. 

When performed using paper maps, the above process is difficult and 
quite time consuming . A database of scanned raster map images, alone, can be 
used to reduce the time and effort expended on step 3 above. However, 
georeferenced raster map images can also be used to eliminate virtually all work 
from steps 2 and 4 above. Since those are generally the most time consuming 
steps under the current process, the value of georeferenced flood maps is 
considerable. 

Using georeferenced flood map raster images, steps 2 and 4 above, are 
replaced by: 

2. A computer system combines the pre-designated outlines of the 
raster map and the georeferencing information to obtain a polygon expressed in 
terms of latitude and longitude that outlines the region included in each flood 
map. Then the system determines which of the polygons contain the address in 
question, which is done using a "point-in-polygon" algorithm. At the conclusion 
of this process, the computer system has identified a map panel (or perhaps a 
small number of map panels) that contains the address. 

4. Since the latitude and longitude of the property are known (by 
virtue of a geocoding phase), the computer system can use the georeferencing of 
the map panels to locate the property on each of the panels found above, thus 
largely eliminating any need for he map analyst to scan the flood map for the 
address location 

Therefore, it is easily seen that by using georeferenced raster maps instead 
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of paper maps, the process of determining flood zones is substantially automated 
and much more efficient. 

It is important to note that while the present invention has been described 
in the context of a fully functional data processing system and/or network, those 
5 skilled in the art will appreciate that the mechanism of the present invention is 

capable of being distributed in the form of a computer usable medixun of 
instructions in a variety of forms, and that the present invention applies equally 
regardless of the particular type of signal bearing medium used to actually carry 
out the distribution. Examples of computer usable mediums include: nonvolatile, 
10 hard-coded type mediums such as read only memories (ROMs) or erasable, 
electrically programmable read only memories (EEPROMs), recordable type 
mediums such as floppy disks, hard disk drives and CD-ROMs, an transmission 
type mediimis such as digital and analog communication links. 

While the invention has been particularly shown and described with 
1 5 reference to a preferred embodiment, it will be understood by those skilled in the 

art that various changes in form and detail may be made therein without departing 
from the spirit and scope of the invention. 
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